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We study multistrange hadrons produced in NA49 and STAR experi- 
ments at center of mass energies varying from y / ijv r /v = 7.61 GeV to 200 
GeV. We show that the yields of S, S and <j> can help to constrain the phys- 
ical conditions present in the hot dense fireball source of these multistrange 
hadrons created in heavy ion collision. We address the question of chemical 
equilibrium of individual quark flavors before and after hadronization and 
offer a few predictions for LHC. 

PACS numbers: 24.10.Pa, 12.38.Mh, 25.75.-q, 13.60.Rj 

1. Introduction 

The Statistical Hadronization Model (SHM) successfully describes par- 
ticle production in heavy ion collisions for a wide range of reaction energies 
studied in the past decade at SPS and RHIC, and has been widely accepted 
as a tool in the analysis of physical phenomena. Special features of multi- 
strange particle yields observed in these experiments inspired this work. 

Particular particle ratios are sensitive to the valence quark content of 
the particles under consideration. For this reason we can explore the hot 
dense fireball flavor content. We gain insight regarding the mechanisms of 
QGP breakup by considering how simple assumptions about dynamics of 
the hadronization alter the result. We extend here our discussion of how 
the yields of multistrange hadrons can be used in such an analysis pQ. 

We begin with a summary of experimental results considered in section 
[5] and introduce the statistical hadronization models in section [3) Our main 
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Fig. 1. Data points of E/<j) Eq. ©, H/K Eq. © and Z/ir Eq. (JTTJ). The straight 
dashed line shows average £/</> — 0.281. 



interest here will be the ratio of 'E(ssq) with 4>(ss), which we address in detail 
in subsection 13.31 This experimental data will help us establish that there 
exists a common set of hadronization conditions for quite different reaction 
conditions, such as centrality, and collision energy. Based on the universality 
of SPS and RHIC results we propose in section 0] several predictions for LHC 
energies. 

2. Experimental results 

Data of interest for us obtained in Pb-Pb collisions by NA49 collabora- 
tion at SPS and in Au-Au collisions by STAR collaboration at RHIC are 
summarized in table [TJ Using these results we obtained the generalized ra- 
tios ofE/(j) (see Eq. ([5|)), E/K (see Eq. @) and E/tt (see Eq. (fTT]) ) shown in 
figure [TJ We note that the ratio of H /(f) is practically constant for all sys- 
tems studied with center of mass energies varying from ^snn = 7.61 to 200 
GeV and number of participants ranging from N part ~ 20 (most peripheral 
collisions) to 350 (most central). 

We are in particular interested in recent data from RHIC obtained as 
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Table 1. Values for the ratios of E/<f> Eq. ©, 3/K Eq. ©, S/tt and Eq. ([IT]) 
obtained from the experimental data of individual particle yields and the resulting 
estimated uncertainty in j q and "f s respectively. Symbol 'E' in the error column 
signals that the particle ratio is the result of the interpolation and/or extrapolation 
needed to account for different centrality bins. 



Experiment 


Ref. 


Centrality 


1OH/0 




10 2 H/K 


hs 


STAR 62 


El SI 


0-5% 


3.04 


E 


4.19 


9.6% 


STAR 62 




5-10% 


3.00 


E 


4.08 


9.2% 


STAR 62 




10-20% 


2.94 


E 


4.06 


9.3% 


STAR 62 


0E1S] 


20-40% 


2.88 


12.5% 


3.79 


E 


blAR 62 


foi pTi rn 
|2J|3J|4] 


40-60% 


2.85 


14.67o 


3.38 


E 


STAR 62 


mmm 


60-80% 


2.49 


19.3% 


2.84 


E 


STAR 200 


[2, 5, 4j 


0-20% 


3.02 


11.8% 


4.06 


12.9% 


SPS 80A 


[amis] 


7% 


3.33 


24.5% 


3.04 


22.7% 


SPS 40A 


IS 13 IS] 


7% 


2.45 


42.1% 


1.89 


18.0% 


SPS 30A 




7% 


2.57 


66.5% 


1.85 


24.3% 


Experiment 


Ref. 


Centrality 


10 3 H/vr 


10 0/K 


lO 2 0/vr 




STAR 62 


01311 


0-5% 


6.22 


1.38 


2.04 




STAR 62 




5-10% 


6.20 


1.36 


2.06 




STAR 62 


mmm 


10-20% 


5.98 


1.38 


2.04 




STAR 62 


mmm 


20-40% 


5.48 


1.32 


1.91 




STAR 62 


mmm 


40-60% 


4.65 


1.18 


1.63 




STAR 62 


mm® 


60-80% 


3.45 


1.14 


1.38 




STAR 200 


[2, 5, 4] 


0-20% 


6.04 


1.34 


to 54+°- 2i 




SPS 80A 


mmm 


7% 


2.6 


0.83 


0.88 




SPS 40A 


mmm 


7% 


3.23 


0.78 


0.83 




SPS 30A 


mmm 


7% 


2.1 


0.63 


0.72 





' For STAR 200 4>/tt considering figure 14 in [2] we give an average of data 
for centrality up to 50%. 



a function of centrality at ^/snn = 62.4 GeV. We performed a centrality 
dependence study shown in figure [2j For every particle species the yield can 
be interpolated by the following form 

/(iV part ) = a-iV p 6 art + c. (1) 

This function is a good empirical approximation of the individual particle 
yields. From the functional form of Eq. ([1]) seen in figure [21 one can see that 
the individual S and 4> yields change by a factor of ten for Au-Au collisions 
at 62 GeV studied in the STAR experiment. When we look at the higher 
energy at STAR and the lower energies at NA49, these individual yields 




vary even more. However, the ratio (as denned by Eq. ([5|)) remains 
within the small range of 0.249 < S/0 < 0.304. 

3. Statistical hadronization considerations 

3.1. Equilibration during hadronization 

Within the SHM approach particle yields are obtained by evaluating the 
accessible phase space. The particle density can be written as, 




where T is the chemical freeze-out temperature, g is the spin degeneracy 
factor, 7 is the phase space occupancy, A is the fugacity factor based on the 
quark content of the given hadron species, and K2 is the modified Bessel 
function of the second kind. 

There are several popular models of statistical hadronization. They 
differ by the assumption made about how the quark content of the QGP 
phase is imaged into the valence quark content of the final hadron state. To 
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begin, one (tacitly) assumes chemical equilibrium of quarks in the quark- 
gluon plasma phase. This equilibrium condition does not have to (and in 
fact it is argued in the following that it cannot) become chemical equilibrium 
of the emerging hadron gas. 

To allow for varying degree of chemical equilibration among final state 
hadrons one introduces the phase space occupancies j q and 7 S which control 
the overall number of quark anti-quark pairs of the light and strange quarks 
respectively. The fugacity factors X q and A s control the relative chemical 
equilibrium (i.e., the relative yields of u,d,s quarks and their respective 
anti-quarks) . 

(i) : In chemical equilibrium, slow hadronization is assumed. During this 
phase transition, all the flavors of quarks have enough time to re-equilibrate 
into their respective hadron gas equilibrium j q = 1 and j s = 1 (q = u,d). 

(ii) : Semi- equilibrium model takes into account that lighter quarks could 
equilibrate faster in hadronization process than strange ones and thus only 
7 S ^= 1, while is fixed to 7 g = 1. 

(iii) : Chemical non- equilibrium assumes a rapid hadronization so that not 
even light quarks have enough time to re-equilibrate. Both strange and 
light quark phase space occupancies are allowed to adopt values different 
from unity as appropriate to describe the observed particle yields. Note 
that when using this most general approach, it is possible to deduce the two 
earlier models as result of data analysis. 

3.2. Multistrange particle production in SHM and model parameters 
To be specific, the directly produced yield of <j) can be expressed as 




It is important to see how fast this series converges. It depends not only 
on the mass to temperature ratio, but also on the phase space occupancies 
and fugacities as they may be greater than unity and appear in increasing 
powers in the series. To get a sense for how fast this sum converges, let us 
examine the following case (table [2]). 

When we calculate the partial contributions to the particle yields of 
interest, we find as shown in table El that the heavier the particle is, the 
faster the sum converges. We have to be careful with less massive particles, 
such as pions, and even kaons. Considering the first seven terms for pions 
gives us only 96% of their yield. For Kaons, we have to take the first 
two terms to get 99.8% of their yield. However, for particles whose mass 
m 3> 100MeV/c 2 , it is sufficient to take the Boltzmann approximation, i.e., 



G 
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Table 2. Example of parameters used to study particle densities. 



Volume 


V 


1000 fm 3 


Temperature 


T 


140 MeV 


u, d Phase space occupancy 


7? 


1.6 


s Phase space occupancy 




2.2 


Fugacities 




1 



just the first term n = 1 in Eq. (J2J: 



N 4-7T , 2 — ( m 



F = s W AmT MrJ' (4) 

One has to be even more careful when using expansion of the Bessel function 
in m/T which converges much slower. 

Table 3. Contribution of individual terms in the sum Eq. ([2]) to the particle yields. 







K~ 


7T 


m[MeV/c 2 ] 


1321 1020 


492 


139 


m/T 


9.4357 7.2857 


3.5143 


0.99286 


N 


0.9862 5.6488 


51.4969 


219.911 




Fractions of contributions order by order: 




n=l 


0.9947 0.9992 


0.9692 


0.6861 


n=2 


1.98 x 10~ 4 1.04 x 10~ 3 


2.89 x 10~ 2 


0.138 


n=3 


6.46 x 10~ 8 1.81 x 10~ 6 


1.52 x 10~ 3 


5.74 x 10~ 2 


n=4 






3.14 x 10~ 2 


n=7 






9.78 x 10~ 3 



3.3. E/4> ratio 

Experimentally measured particle yield consists not only of directly pro- 
duced particles, but also a contribution from heavier resonances decaying 
into it and thus increasing its apparent yield. In general the experimental 
results comprise these contributions. Therefore, we include in our calcula- 
tions directly produced H's and </>'s and the contributions from decays of 
heavier resonances. The most significant contribution comes from H*(1530) 
with spin degeneracy g = 4. We use the decay tables as implemented in the 
SHARE program [TT)] . 

Particle yield ratios are very sensitive to the valence quark content and 
we can take advantage of this by considering specific ratios. One uses ratios 
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since when one takes the same number of hadrons in the numerator as in 
the denominator, the overall normalization is canceled. We also consider the 
product of a particle and its antiparticle in order to eliminate the chemical 
potentials, which are opposite, and therefore the fugacity factors also cancel. 
Finally, by the choice of hadrons involved, one can change how the ratio 
depends on the phase space occupancies. 

For these reasons, we look first at the generalized ratio of 3 to </>, defined 

as: 



The ratio of is proportional to the light quark phase space occupancy 
7 g and a function / of temperature T, 



which in the Boltzmann approximation is determined by degeneracies gi 
and masses of the hadrons contributing to the yields and over which we 
sum using the Boltzmann approximation. 

As evident from Eq. ([2]), the effect of quantum statistics is decreasing 
with the hadron mass m when compared to the freeze-out temperature 
T. We calculated the magnitude of this effect to be 0.25% for the S/0 
ratio and is therefore negligible. One has to remember that the formula for 
the residual function Eq. (0) does not account for strong decays in which 
flavor content is redistributed, nor does reflect any possible weak decays 
influencing the yields in question. 

Global fits [111 [12"| [T3] of the systems under consideration show that all 
these systems hadronize at the same temperature T in the non-equilibrium 
model. This, together with the observed constant ratio of implies com- 
mon hadronization conditions for all the systems studied. Another immedi- 
ate implication is that the light quark phase space occupancy 7 9 is constant 
and has the same value for all the energies, centralities and systems studied. 

To obtain the results presented here we actually use the SHARE program 
[lOj . We consider several values of the S/(/> ratio and plot the results in the 
7q-T plane in figure [3l In combination with the global fit results, there 
are two possible regions to look for possible hadronization conditions. The 
first region corresponds to light quarks being in chemical equilibrium. This 
is represented by 7 9 = 1, and from figure O we see temperature T ~ 170 
MeV. To locate the second region, we allow for fast hadronization and relax 
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Fig. 3. Lines of a constant ratio 10S/</> £ [0.8,4.5] (as defined in Eq. in the 
T-7 9 plane. The lines for j q = 1 and 7 9 = 7™' are presented by solid black lines. 
The thick gray line represents the average result of all SPS and RHIC. 



light quark phase space occupancy j q 7^ 1 . When varying the magnitude of 
jq, one has to remember that there is a critical value of this parameter for 
which pions begin to form a Bose-Einstein condensate. This critical value 
is given by the mass of the lightest, non-strange meson (tt°) by: 



7 9 <l q 



crit 



exp 



(7) 



This is shown as an upper limit in figure [3j 

In conclusion of this analysis, we observe that it is possible for different 
reaction systems to occupy different values of T and 7^ along the solid black 
line in figure [3l However, much more plausible interpretation is that all 
reaction systems considered hadronize at the same set of T,j q . This result 
is strongly supported by the global fit of all hadron yields we performed and 
in which near critical value of 7^ ~ 1.6 is favored. Figure then implies 
temperature T ~ 140 MeV. 
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3.4- Hadronization of strangeness 

In the previous subsection, we found that the ratio S/0 is proportional 
to light quark phase space occupancy 7 9 . To study strangeness occupancy 
7s we have to identify a ratio proportional to 7 S . By the same reasoning as 
before, we are motivated to examine the ratio of S to K defined as: 

It is evident from both tableland figure [1] that S/K is not constant. 
The Au-Au, y/s^N = 62.4 GeV, RHIC data which inspired this study show 
a systematic variation with a factor of 1.5 difference between central and 
peripheral collisions. When we consider the other systems and energies 
studied, the factor increases to more than two. We believe this is due 
entirely to variation in 7 S since as we argued that T = Const. 

When 7 S is allowed to vary, we have to remember, especially in the 
view of the expected conditions at LHC where 7 S could turn out to be 
large, that there is also a critical value of 7 S . It turns out that it is the rj 
meson (rj = 0.45(uu + dd) + 0.55ss) which is the first to condensate when 
7s increases. The critical value for 7 S is a function of temperature and can 
be obtained from the condition: 

0.45 7 2 + 0.557s 2 = exp (^j . (9) 

The critical value of "y s is not as sensitive to 7 g as it is to temperature. For 
the two scenarios discussed, T = 140 MeV and T = 170 MeV, we get the 
critical value of 7s"* = 9.47 and 7s"* = 6.68 respectively. 

4. Predictions 

4-1. General considerations 

We have argued in the last section, that the temperature is constant 
for all the systems, and we have attributed the change in relative yield of 
E and K to the change in j s . This implies that the strangeness initially 
present in the fireball of QGP does not have enough time to re-equilibrate 
during the hadronization, and the strange phase space in the hadron gas 
phase can be under- or over-populated. As a consequence, we conclude that 
chemical equilibrium model (i) in subsection 13.11 which assumed both light 
and strange quarks to be in chemical equilibrium after hadronization, is no 
longer a viable description of particle yields in heavy-ion collisions. 
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With the range of 7<j,7 s and T established by measuring 3/0 and S/K, 
we can now predict the ratios of particles which are harder to observe ex- 
perimentally. According to the quark content of the particles in question, 
we establish how their ratio depends on phase space occupancies and with 
that predict the behavior and possible values of their ratio. 

Both the ratios S/K and f2/0 are proportional to 7 S . In the top section 
of figure H] we show fixed values of 3 /K and the band of experimental data 
for this ratio in the T — j s plane. In the bottom frame of figure H] we show 
the corresponding f2/0 ratio and find a narrow range: 



Similarly, with temperature T, 7 g and j s constrained to a small range 
of values, we can predict other ratios of particles, such as (ft/tr or 3/n (see 
figure [Q) defined by 



which are both very sensitive to the amount of strangeness in the emerging 
hadron gas due to their dependence on 7^. 



One very simple and solid prediction we can make based on the univer- 
sality of the results for 3/0 is that 3/0 = 0.28 ± 10% at LHC too. This 
universal hadronization condition for T, 7^ has been observed so far in both 
SPS and RHIC results. We strongly believe that this yield will be observed 
in heavy ion reactions at LHC. 

We also expect considerably greater absolute yield of at LHC, since 
one can anticipate that there is a greater strange quark pair density produc- 
tion because of more extreme initial conditions [13] • We found in a kinetic 
exploration of strangeness production that the strangeness over entropy ra- 
tio at LHC is at least s/S — 0.037. This implies, as shown in [15J . that the 
ratio of Js/^fq — 1-55. We use this value in our further LHC predictions. 
Here we note that since the size of reaction volume changes, as does baryon 
and strangeness content, several SHM parameters such as V, 7 s ,A g ,A s are 
expected to change. Thus the hadronization universality only refers to T, 7^. 

We show the required value of 7 S as a function of the measured value of 
4>/tt in figure[S]for both chemical semi-, and non-equilibrium. Also we show 
(colored) bands of data from different experiments and also a prediction 




;+ 



(10) 





4-2. LHC conditions 
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Fig. 4. Lines of constant ratio S/K (top) and Q,/<f> (bottom). Experimental data 
from most central 0-20% STAR 62 are indicated by a thick line (top), and is 
assumed in the bottom frame as a prediction. See text for more detail. 
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Fig. 5. 7 S as a function of the relative tp/n yield Eg. (jlip in two hadronization 
scenarios, see text. The horizontal solid black line shows the chemical equilibrium 
with 7 S = 1. For experimental data see tabled] 

for LHC results. It is evident from the figure that the temperature does 
not have a great effect on this ratio when we fix 7^ = 1, a very fortuitous 
result — the increase of 0-yield with T is nearly compensated by increase 
in 7r-yield. Thus, if at LHC cfi/ir > 0.03 we can be sure that irrespective of 
the model of hadronization the value of j s > 1. 

7 S > 1 cannot be build up from a small initial value in hadron based 
strangeness production since chemical equilibrium is reached at 7 S = 1 . 
On the other hand, when the source of hadrons is a strangeness rich initial 
phase, such as chemically equilibrated quark-gluon plasma, hadronization 
can lead to phase space overpopulation in final state. This is true in par- 
ticular for a fast hadronization process, in which strangeness content does 
not have time to re-annihilate. 

5. Conclusions 

The experimental results from SPS and RHIC provided an important 
observation, that the ratio S/<^ is nearly constant for wide range of energies, 
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centralities and even across different systems, while the individual yields 
change by more than an order of magnitude. As this ratio is governed 
by light quark phase space occupancy j q , and hadronization temperature 
T, than these conditions are very likely universal for all the experiments. 
This universality in regard to T, ^ q is strongly supported by global data 

fits na ng. 

A different ratio, H/K, which is proportional to strange phase space 
occupancy 7 S , changes by a factor of 2, which implies that 7^ changes sig- 
nificantly between different experiments and centralities. Thus, we learn 
that the full chemical equilibrium hadronization is not a viable model. 

We further find that some experimentally observed yields of multistrange 
hadrons, such as cp, are only compatible with 7 S > 1. For this to hap- 
pen, we require strangeness rich initial phase (quark-gluon plasma) and 
fast hadronization, so that strangeness abundance does not have enough 
time to re-equilibrate during hadronization. 

Our study allowed us to predict certain particle ratios for which precise 
data is not available at present experiments and for the upcoming LHC. For 
example, we expect 0/0 G (5.5 x 1(T 2 ,7.0 x 1(T 2 ) at RHIC 62 GeV. We 
expect the ratio 3/(f> to remain in the same range of values at LHC as it 
has been measured at SPS and RHIC, i.e., E/<p G (0.249,0.304). 
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